High Angular Momentum Core States and Impurity Effects in the Mixed State of 

d-Wave Superconductors 
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The local quasiparticle spectra around the vortex core of a d-wave superconductor in the mixed 
state is studied by solving the lattice Bogoliubov-de Gennes equations self-consistently. It is shown 
that in addition to the zero-energy states, there also exist core states of high angular momentum. 
A nonmagnetic unitary impurity sitting at the core center is nondestructive to these core states and 
the zero-energy resonant state induced by itself is still visible in the local quasiparticle spectrum. 
The calculated imaging shows a fourfold "star" shape of the local density of states whose orientation 
is energy dependent. It is also found that, although the zero-energy states are extended, the core 
states of finite energy level could be localized. 
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In the Abrikosov vortex state (also called mixed state 
or Schubnikov phase) of type-II superconductors, the 
magnetic flux penetrates into the system in quantized 
units $o = hc/2e. In the region of vortex core, the su- 
perconducting order parameter (or pair potential) is con- 
siderably suppressed from its zero-field value A n , going 
to zero at the vortex core center at a distance of several 
coherence lengths £o ■ This reduction of the order parame- 
ter can affect various physical properties of the supercon- 
ductor, for example, the nature of local electronic states 
around the Fermi energy. The vortex core problem of 
an s-wave superconductor was first studied in the classic 
papers of Caroli, de Gennes, and Matricon 0]. By solv- 
ing approximately the Bogoliubov-de Gennes equations 
for an isolated vortex line in the limit of k = A/£o 3> lj 
where A is the magnetic penetration depth, these au- 
thors showed the existence of the low-lying bound quasi- 
particle states inside an s-wave vortex core. Later on, 
Bardeen et al. Q extended the calculation to all values 
of k by determining the pairing potential and magnetic 
field with a variational expression for the free energy, and 
the qualitative conclusion is similar to that of Caroli et 
al. Recent theoretical attempts ]|, |[ |5[ [| 0, |) at un- 
derstanding the electronic states of vortex lines in type- 
II superconductors was revived by STM experiments on 
NbSe2 |9) , which not only showed a strong enhancement 
of the zero-bias tunneling conductance as a manifesta- 
tion of the presence of bound states in the vortex core, 
but also provided images of vortex. For an s-wave su- 
perconductor, the energy gap open at the Fermi surface 
is a constant. One can naively think that, the spatial 
variation of the pair potential is analogous to a poten- 
tial well for the quasiparticle, with depth Ao and radius 
£0 ■ The energy levels of the quasiparticle bound states is 
roughly the multiples of Eq = fi 2 /2m£ 2 w A^/Ep, where 
the energy is measured relative to the Fermi energy Ep. 
For a d-wave superconductor as recently established in 
Cu02-based high-T c cuprates, the energy gap has a max- 
imum along the copper oxide bonds and is closed along 



the diagonals of square Cu02 lattice. In some respect, 
the spatial variation of the pair potential is analogous 
to a potential well carved with four notches along the 
nodal directions, so that quasiparticles can leak through 
the notches. It was shown for the first time by Wang 
and MacDonald [[To) that there appears a single broad 
peak at zero energy in the local density of states (LDOS) 
at the center of a pure d-wave vortex core so that the 
quasiparticle states are virtually bound. Recently, the 
split-peak structure around zero bias in the local differ- 
ential tunneling conductance at the vortex core center in 
YBaaCuaO^ @ and Bi 2 Sr 2 CaCu 2 8+5 @, [l§ has 
been observed by scanning tunneling microscopy (STM) . 
These experimental observations stimulated intensified 
study of the quasiparticles in the vortex core of high- 
T c cuprates @, [l§ £§ [l§ ^ ^ |§, || || g f§, 
leading to a variety of scenarios for the explanation of 
experimental data. No consensus on the mechanism has 
emerged. However, on the other hand, additional fea- 
tures of quasiparticle states of a pure (i-wave vortex have 
not been well addressed even in the framework of BCS- 
type description. Specifically, nearly all of the theoreti- 
cal and experimental works on the quasiparticle states of 
d-wave vortices focus on the LDOS directly on the core 
center, whether there exists high angular momentum core 
states have not been paid much attention. Furthermore, 
the imaging of the tunneling conductance of a d-wave 
vortex, as in the impurity case |2(| , has not been experi- 
mentally available. In this paper, we present a complete 
study of the electronic structure around one of the d-wave 
vortices forming the Abrikosov square lattice by solving 
self-consistently the Bogoliubov-de Gennes (BdG) equa- 
tions within a tight-binding model, and new features are 
exposed for the first time. 



The electronic structure of a d-wave superconduc- 
tor can be described by the quasiparticle wavefunctions 
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which satisfy the BdG equations 
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Here the summation is over the nearest neighbor 
sites. The single particle Hamiltonian reads 7i%j = 
—te llpij 5i +T ,j + (Ui— n)5ij where t is the hopping integral, 
fj, is the chemical potential, f/j is the single particle poten- 
tial describing, if any, the effects of impurities, defects, or 
crystal field, and r = ±e x , ±e y are the unit vectors along 
the crystalline x and y axes, respectively. In the mixed 
state, the magnetic field effect is included through the 
Peierls phase factor ip^ = j r % A(r) • dr. By assuming 
the superconductor under consideration is in the extreme 
type-II limit where k goes to infinity so that the screen- 
ing effect from the supercurrent is negligible. The vector 
potential A can then be approximated by the solution 
V x A = Hz where H is the magnetic field externally 
applied along the c axis. The enclosed flux density within 
each plaquette is given by ^ n (fij — ^g" . Notice that 
the quasiparticle energy is measured with respect to the 
Fermi energy. is the spin-singlet d-wave bond pair 
potential, which is subject to the self-consistency condi- 
tion: 
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with V the strength of nearest neighbor effective attrac- 
tion between electrons. Hereafter we measure the length 
in units of the lattice constant a and the energy in units 
of the hopping integral t. Within the Landau gauge the 
vector potential can be written as A = {—Hy, 0, 0) where 
y is the y-component of the position vector r. We intro- 
duce the magnetic translation operator T mn r = r + R 
where the translation vector R = mN x e x + nN y e y with 
N x and N y the linear dimension of the unit cell of the 
vortex lattice. To ensure different T mn to be commutablc 
with each other, we have to take the strength of mag- 
netic field so that the flux enclosed by each unit cell has 
a single-particle flux quantum, i.e, 2<& . Therefore, the 
translation property of the superconducting order param- 
eter is A(T mn r) = e lx ( r ' R ) A(r) where the phase accu- 
mulated by the order parameter upon the translation is 
x(r,R) = ^A(R)-r— Amnn. From this property, we can 
obtain the magnetic Bloch theorem for the wavefunction 
of the BdG equations: 
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Here f is the position vector defined within a given 



unit cell and k 
0, 1, . . . , M x>y 



M^N^ X ^ M y N y ^V with l*,V 

1 are the wave vectors defined in the first 
Brillouin zone of the vortex lattice and M X N X and M y N y 
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FIG. 1: The local density of states as a function of energy 
at various distance away from the d-wave vortex core center 
(0, 0). The coordinate is measured in units of the lattice con- 
stant. The solid line in each panel shows the result in the 
clean limit. The result with a nonmagnetic unitary impurity 
sitting at the core center are shown by the dashed line. The 
bulk density of states in the absence of the magnetic field and 
impurities is also plotted for a reference by the dotted line. 
For the chosen unit cell size 42 x 21, the site (10,0) is the 
midpoint between two nearest neighboring vortices along the 
x direction. The other parameter values: V = 1.2, nj — 0.85, 
and T = 0.01. 



are the linear dimension of the whole system. We use ex- 
act diagonalization method to solve the BdG equation (|l|) 
self-consistently. The thermally broadened local density 
of states is then evaluated according to 
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where a prefactor comes from the spin degeneracy, and 
f'(E) is the derivative of the Fermi distribution function 
f(E) — l/[exp(E/k B T) + 1] with respect to the energy 
E, N c = M x x My is the number of magnetic unit cells. 
Pi(E) is proportional to the local differential tunneling 
conductance at low temperatures which could be mea- 
sured by STM experiments j27|. 

In the following calculation, we take V = 1.2, and 
the filling factor 0.85 by adjusting the chemical poten- 
tial. The temperature is chosen to be T = 0.01. With- 
out loss of generality, we consider the unit cell of size 
N x x N y = 42 x 21, and the number of the unit cells 
N c = 21 x 42. This choice will give us a square vortex lat- 
tice. The numerical calculation shows as expected that 
each unit cell accommodates two superconducting vor- 
tices each carrying a flux quantum $o- The spatial vari- 
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ation of the d-wave order parameter is similar to Fig. 1 
of Ref. 1 10 : It decreases continuously to zero from the 



its bulk value as the vortex core center is approached 
in the length scale £ — ^ w _F/2A maa . rs 5, with the de- 
pleted region extending farther in the diagonal direction 
of the square lattice. Here A max ss 0.38 is the gap edge 
defined as the energy position at which the bulk DOS 
becomes maxima, and vf is the Fermi velocity. For the 
given model parameter, the induced extended s-wave or- 
der parameter around the vortex is so small that its effect 
on the LDOS should be negligible. In Fig. |l| we plot the 
LDOS as a function of energy at various distance from 
the vortex core in a d-wave superconductor. The solid 
line in each panel represents the result of a clean sys- 
tem, while the dashed line for the case that a unitary 
impurity (Uq — 100) sits at the vortex core center. For 
comparison, we have also displayed the bulk LDOS (the 
dotted lines) in the absence of magnetic field and im- 
purities. We choose the core center as the coordinate 
origin (0,0). For the given unit cell size, the site (10,0) 
is the midpoint of two nearest neighboring vortices along 
the x direction. Notice that the three curves in the right- 
bottom panel almost coincide with each other, indicating 
that the LDOS spectrum at the midpoint has recovered 
the bulk DOS. The asymmetry of the LDOS spectrum 
with respect to the Fermi energy comes from the break- 
ing of the particle-hole symmetry. In the absence of the 
unitary impurity, the LDOS at the core center i = (0, 0) 
shows a single resonant peak around the Fermi energy. 
Although the spectrum at the center does not exhibit 
a split-peak structure, the suppression of the coherent 
peak at the gap edge A max seems to be consistent with 
the experimental observations. At this point, we also 
calculate the LDOS at the core center with various val- 
ues of the bulk order parameter and find that the broad 
zero-energy peak always shows up, which leads us to con- 
clude that the quasiparticle excitations along the nodal 
direction on the Fermi surface should be the origin. At a 
finite distance away from the core center, the single peak 
in the LDOS at the core center continuously evolves into 
two peaks, which are located at ±5E(i). The spacing 
between these two peaks, i.e., 2SE(i) increases with the 
distance away from the core center. In addition, the am- 
plitude of these peaks decreases with the distance from 
the core center and finally they merge into the contin- 
uum of scattering state with E > E max . We remark 
that, since the induced extended s-wave order parame- 
ter is almost negligible, this origin for the obtained split 
peaks should be excluded. Considering the distance de- 
pendence of these subgap peaks in the LDOS, we believe 
instead that they result from the core states with a high 
angular momentum. The effect of a nonmagnetic impu- 
rity in the unitary limit is twofold: First, it plays the 
role of a pinning center. We introduce such an impu- 
rity into the unit cell randomly. It is found that when 
this impurity has a distance from the vortex core center 



within the range of a coherence length, the vortex will 
be dragged onto the impurity site. The order parameter 
profile of the pinned vortex is quite similar to that of a 
clean case. Out of this range, the pinning effect is rather 
weak. Second, it is both theoretically ^8) and experi- 
mentally (2^, Q well established that a nonmagnetic 
impurity itself will induce resonant states in a d-wave 
superconductor in contrast to its counterpart in a con- 
ventional s-wave superconductor. For the convenience 
of discussion, we consider the lattice as a bipartite, i.e, 
those sites i = (ix,iy) with index having even (— \y x + l v 
form sublattice A, which contains the core center, while 
those with index having odd (— l) KB +*i' form sublattice B. 
The panels in the left column of Fig. |l| display the LDOS 
at the sites belonging to A, and the panels in the right 
column show the LDOS at the sites belong to B. When 
such an impurity pins the vortex, the LDOS spectrum at 
the core center vanishes due to the strong potential scat- 
tering from the impurity [See the dashed line as shown 
in the panel of Fig. |l| corresponding to the site (0,0)]. 
However, the intensity of the single zero-energy peak in 
the LDOS at the nearest neighbor site of the core cen- 
ter is strongly enhanced. This enhancement comes from 
the additional contribution of the near-zero-energy res- 
onant state induced directly by the impurity itself p8| . 
As the measured point is away from the core center, the 
LDOS exhibits no central peak (near the Fermi energy) 
at the sites belonging to the sublattice A. In particular, 
the density of states at zero energy is more depressed on 
these sites as compared to the impurity free case. Nev- 
ertheless, the LDOS on the sites belonging to the sub- 
lattice B has a zero-energy peak, the intensity of which 
decays with the increased distance away from the core 
center. This pattern of the LDOS around the zero en- 
ergy is solely caused by the impurity-induced resonant 
state, which has an oscillation behavior 0, |TJ. Also in- 
terestingly, it can be seen that the location of the LDOS 
peaks associated with the high angular momentum core 
states is quite insensitive to the impurity scattering. In 
this sense, we can say that the unitary impurity is non- 
destructive to the vortex from the point of view of core 
states. The above effect of the single nonmagnetic uni- 
tary impurity on the electronic structure of a d-wave vor- 
tex is dramatically different from its effect in the case of 
an s-wave vortex. Since there is no resonant states in- 
duced around this impurity in an s-wave superconductor, 
it drives the energy position of the low-level core states 
away from the Fermi energy, and simultaneously the in- 
tensity of the corresponding peaks in the LDOS becomes 
weaker with the impurity strength p2| |. 

To take a closer inspection of the nature of the elec- 
tronic structure around a d-wave vortex, we calculate the 
spatial distribution of the LDOS around the vortex core 
at the fixed energy E\ = and E2 = 0.15. Si is the peak 
position of the LDOS at the core center, E2 is that at the 
site (2,1). For the clean case, we find strong anisotropy of 
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the spatial distribution of the LDOS, which has a fourfold 
"star" shape. At zero energy (E\ = 0), the magnitude of 
the LDOS has tails extending along the diagonals of the 
square lattice, and decays rapidly along the Cu0 2 bonds. 
This behavior can be understood as follows: When the 
energy gap vanishes along the nodal directions, the qua- 
siclassical orbital of these low-lying quasiparticle states 
is not closed. Therefore, the zero-energy states are ex- 
tended along these directions and consequently, in con- 
trast to the s-wave case, they are not truly bound core 
states. However, at a finite energy (E 2 = 0.15), the situ- 
ation is reversed and the LDOS decays more slowly along 
the Cu02 bonds than along the diagonals of the square 
lattice, that is, the "star" shape is rotated by almost 
7r/4. In particular, the decay even along the bonds is so 
short ranged that the core states at this energy is approx- 
imately localized. This energy-dependent orientation of 
the "star" shape around a d-wave vortex arises from the 
d-wave pairing symmetry itself, which is in striking con- 
trast to the s-wave case, the shape rotation was explained 
by introducing a perturbation term ^ . When the impu- 
rity is added on the core center, the above conclusion is 
not changed except that the LDOS directly at the impu- 
rity site and the "star" becomes more extended. 

To summarize, we have presented a study of the quasi- 
particle states around the vortex core of a d-wave super- 
conductor in the mixed state. We have shown that, in 
addition to the zero-energy core states, there also exist 
core states with high angular momentum. It has also 
been found that a nonmagnetic impurity is not destruc- 
tive to these core states, and the impurity-induced near- 
zero energy resonant state [£8| can still be visible in the 
LDOS. The calculated images demonstrate that all core 
states are not truly localized. So far, it seems that all 
theoretical and experimental studies focused on the local 
electronic structure at the vortex core center, on which 
much controversy just arises. The test of the existence 
of the high angular momentum core states, as well as the 
effects of the single atomic impurity, is now accessible to 
the experiments. Currently, the STM is the technique 
of choice that allows the observation of vortex states. It 
has recently been improved in both the spatial (atomic 
scale) and energy (sub-meV) resolution |26|. We suggest 
that the STM tip be scanned horizontally in the region 
of several lattice constants around the vortex core center 
to observe the predicted changes in intensity and energy 
position of peaks. If the above predictions are indeed 
experimentally observed, it will help us to understand 
the delicate nature of the electronic structure of high-T c 
vortices. 
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FIG. 2: The spatial distribution of the LDOS around a d-wave 
vortex at the energy Ei = (top row) and E2 = 0.15 (bottom 
row). The left column shows the results in the clean case. 
Those with a nonmagnetic unitary impurity sitting at the 
core center are displayed on the right column. The measured 
size is 21 x 21. The other parameter values are the same as 
Fig. §. 



[2 

[3. 
[1 

[5 

[6 

l? 

[8 

[9 

[10; 

[11 
[12 

[13 
[14 

[is; 

[16 
[17 



C. Caroli, P. G. de Gennes, and J. Matricon, Phys. Lett. 
9, 307 (1964); C. Caroli and J. Matricon, Phys. Kondens. 
Mater. 3, 380 (1995). 

J. Bardeen et al, Phys. Rev. 187, 556 (1969). 
J.D. Shore et al, Phys. Rev. Lett. 62, 3089 (1989). 
U. Klein, Phys. Rev. B 40, 6601 (1989); Phys. Rev. B 
41, 4819 (1990) 

F. Gygi and M. Schluter, Phys. Rev. Lett. 65, 1820 
(1990); Phys. Rev. B 41, 822 (1990); 43, 7609 (1991). 
S. Ullah and A. Dorsey and L. J. Buchholtz, Phys. Rev. 
B 42, 9550 (1990). 

B. Pottinger and U. Klein, Phys. Rev. Lett. 70, 
(1993). 

Y.-D. Zhu, F.C. Zhang, and M. Sigrist, Phys. Rev. 
1105 (1995). 

H. F. Hess et al, Phys. Rev. Lett. 62, 214 (1989); 
Hess et al, ibid. 64, 2711 (1990). 

Y. Wang and A. H. MacDonald, Phys. Rev. B 52, R3876 
(1995). 

I. Maggio-Aprile et al, Phys. Rev. Lett. 75, 2754 (1995). 
S. H. Pan et al, Phys. Rev . Lett. 85, 1536 (20 00). 



2806 



B 51, 



H. F. 



:ond- mat/010552^ . 



B. W. Hoogenboom at al. 

A. Himeda et al, J. Phys. Soc. Jpn. 66, 3367 (1997). 
Y. Morita, M. Kohmoto, and K. Maki, Phys. Rev. Lett. 
78, 4841 (1997). 

D. P. Arovas et al, Phys. Rev. Lett. 79, 2871 (1997). 
M. Franz and Z. Tesanovic, Phys. Rev. Lett. 80, 4763 



5 



(1998). 

[18] M. Takigawa, M. Ichioka, and K. Machida, Phys. Rev. [25] 

Lett. 83, 3057 (1999). [26] 

[19] K. Yasui and T. Kita, Phys. Rev. Lett. 83, 4168 (1999). [27] 
[20] J. H. Han and D.-H. Lee, Phys. Rev. Lett. 85, 1100 

(2000); Q.-H. Wang, J. H. Han, and D.-H. Lee, Phys. [28] 
Rev. Lett. 87, 167004 (2001). 

[21] B. M. Andersen, H. Bruus, and P. Hedegard, Phys. Rev. [29] 

B 61, 6298 (2000). [30] 

[22] C. Wu, T. Xiang, and Z.-B. Su, Phys. Rev. B 62, 14427 [31] 
(2000). 

J. Kishine, P. A. Lee, and X.-G. Wen, Phys. Rev. Lett. [32] 



[23] 

[24] Jian-Xin Zhu and C. S. Ting, Phys. Rev. Lett. 87, 147002 



86, 5365 (2001); fond- mat/0108478 
ridl 



(2001). 

C. Berthod and B. Giovannini, sond-mat /0 1 05425] . 
S. H. Pan at al, Nature 403, 746 (2000). 
M. Tinkham, Introduction to Superconductivity (McGraw 
Hill, New York, 1996). 

M. I. Salkola, A. V. Balatsky, and D. J. Scalapino, Phys. 

Rev. Lett. 77, 1841 (1996). 

E. W. Hudson et al, Science 285, 88 (1999). 

A. Yazdani et al, Phys. Rev. Lett. 83, 176 (1999). 

J.-X. Zhu, C. S. Ting, and Chia-Ren Hu, Phys. Rev. B 

62, 6027 (2000). 

Q. Han, L.-y. Zhang, and Z. D. Wang, Phys. Rev. B 62, 
5936 (2000). 



